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PREFACE 


In  this  Memorandum  the  dependence  of  plane  array  gain  upon 
thinning,  steering  and  element  directivity  is  discussed  and  a 
systematic  way  to  avoid  grating  lobes  in  highly  thinned  arrays  is 
introduced.  The  ideas  are  illustrated  by  some  detailed  examples. 
The  Memorandum  should  interest  people  designing  thinned,  elec¬ 
tronically  steered,  plane  arrays  for  modem  radars. 
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SUMMARY 

Two  related  topics  In  the  design  of  plane  phased  arrays  are 
discussed.  One  is  the  dependence  of  gain  upon  thinning,  steering 
and  element  directivity;  the  other  is  a  systematic  procedure  for 
locating  the  elements  of  highly  thinned,  plane  arrays. 

A  few  general  observations  on  the  effects  of  thinning  and 
steering  can  be  made  from  the  expression  for  the  gain  of  an  array 
of  elements  having  an  arbitrary  directivity  factor.  The  speciali¬ 
zations  to  isotropic  elements  and  to  infinitesimal  dipoles  are  then 
made  and  these  observations  are  verified  for  some  actual  arrays. 

A  systematic  way  to  randomize  element  locations  in  a  thinned 
array  while  still  approximating  some  basic  excitation  taper  is 
also  suggested  and  illustrated.  The  thinning  scheme  is  based  on  an 
equal  area  approximation  to  a  continuous  excitation  distribution, 
with  the  equal  areas  being  laid  out  in  the  form  of  an  expanding 
spiral.  This  procedure  insures  against  accidental  element  align¬ 
ment  in  any  particular  direction.  Its  effectiveness  is  shown  by 
applying  it  to  a  100- element  array  and  plotting  several  sections 
of  the  resulting  directivity  pattern. 
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I.  INTRODUCTION 

Two  related  topic''  in  the  design  of  plane  phased  arrays  are 
discussed.  One  is  the  dependence  of  gain  upon  thinning,  steering 
and  element  directivity;  the  other  is  a  systematic  procedure  for 
locating  the  elements  of  highly  thinned,  plane  arrays. 

A  few  general  observations  on  the  effects  of  thinning  and 
steering  can  be  made  from  the  expression  for  the  gain  of  an  array  of 
elements  having  an  arbitrary  directivity  factor.  The  specializations 
to  isotropic  elements'  and  to  infinitesimal  dipoles  are  then  made  and 
these  observations  are  verified  for  some  actual  arrays. 

A  systematic  way  to  randomize  element  locations  in  a  thinned 
array  while  still  approximating  some  basic  excitation  taper  is  also 
suggested  and  illustrated.  The  thinning  scheme  is  based  on  an  equal 
area  approximation  to  a  continuous  excitation  distribution,  with  the 
equal  areas  being  laid  out  in  the  form  of  an  expanding  spiral.  This 
procedure  insures  against  accidental  element  alignment  in  any  par¬ 
ticular  direction. 

The  dependence  of  gain  on  element  spacing  has  been  described  by 

Milazzo  and  D1 Angelo ^ for  the  case  of  a  uniformly  spaced,  broadside 

(2) 

linear  array.  Bickmore  considers  a  square  current  sheet,  linearly 
polarized  parallel  to  one  side  of  the  square.  For  steering  around 
the  axis  of  polarization  he  derives  the  dependence  of  3  db  beamwidth 
on  steering  and  shows  that  the  cosine  approximation  for  effective 
aperture  holds  to  within  a  beamwidth  or  so  of  endfire.  Von  Aulock's 
expository  article  contains  similar  material  on  the  beamwidths  for 


§ 
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electronically  steered  rectangular  arrays  with  unifoxm  spacing. 
Although  results  like  those  in  Ref.  (l)  are  easily  demonstrated  for 
plane  arrays  also,  the  facts  are  apparently  not  universally  appreci¬ 
ated.  One  difficulty  is  the  number  of  parameters  whose  effects  are 
to  be  disentangled  and  separately  accounted  for.  In  this  Memorandum 
the  effects  will  be  discussed  in  connection  with  the  detailed  example 
which  occupies  about  half  the  pages.  Although  some  qualitative 
statements  can  be  made  from  the  general  expression  for  gain,  the 
actual  fluctuations  in  gain  as  a  thinned  array  is  steered  will  depend 
strongly  on  the  array  configuration  itself.  For  instance,  a  square 
array  with  uniform  spacing  of  a  wavelength  or  so  would  have  very  wide 
gain  fluctuations,  while  the  example  to  be  considered  has  relatively 
mild  gain  fluctuations  with  steering. 

The  example  also  furnishes  a  convenient  way  to  introduce  a  new 

approach  to  the  problem  of  spurious  main  lobe  reduction  in  highly 

(U) 

thinned  plane  arrays.  The  approach  of  King,  Packard  and  Thomas 
is  essentially  intelligent  cut-and-try.  Baklanov,  Pokrovski,  and 
Surdutovich^  apply  a  Newton-Raphson  iterative  procedure  to  solve 
for  the  locations  of  the  N  elements  of  a  symmetrical  linear  array 

(6) 

whose  first  [N/2]  side  lobes  have  a  pre-assigned  level.  Andreasen 
has  also  used  an  iterative  procedure  to  deduce  element  distributions 
for  which  the  nearest  grating  lobes  or  other  large  lobes  are  held 
down  to  a  level  appropriate  to  the  number  of  elements .  All  three 
papers  are  concerned  with  linear  arrays.  For  plane  arrays  or  arbi¬ 
trary  volumetric  distributions  there  have  been  some  statistical 

(7) 

analyses  of  large  random  sets  of  radiators.  Cover's  Memorandum 
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is  typical,  though  not  directly  concerned  with  spurious  lobe  re- 

/Q\ 

duction.  Maher  and  Cheng'  ;  consider  the  statistical  effects  of  re¬ 
moving  a  relatively  small  number  of  elements  from  a  linear  array. 

The  only  reported  explicit  synthesis  of  a  thinned,  plane  array  seems 

{9) 

to  be  that  of  Willey.  The  method  used  in  this  Memorandum  is  very 
close  to  Willey's,  which  is  based  on  the  idea  of  approximating  a 
continuous  aperture  distribution  with  a  discrete  array  of  equal 
radiators  by  partitioning  the  continuous  distribution  into  areas  of 
equal  weight  and  assigning  a  radiator  to  each  such  area.^^ 

It  differs  from  Willey's  in  providing  a  nearly  automatic  assign¬ 
ment  of  element  locations  without  the  need  to  make  individual  choices 
by  hand  from  a  small  number  of  possibilities. 
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II.  GAIN  OF  PLANE  PHASED  ARRAYS 

Let  N  radiators  be  located  at  points,  P^,  in  the  equatorial 
plane,  8  =  n/2,  of  a  spherical  coordinate  system  (see  Fig.  l).  The 
locations  will  be  specified  in  either  rectangular  coordinates, 
(x^,y^),  or  polar  coordinates,  as  convenient.  Let  the 


Fig.  1  —  Coordinates 
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amplitude  and  phase  of  the  excitation  of  the  k^1  element  be  a^  and 
bk,  so  that  the  complex  excitation  is  given  by  or^  *  a^  exp(ib^). 

Consideration  will  be  restricted  to  cases  deacribable  by  a 
scalar  or  one-dimensional  analysis.  These  include  acoustic  radiators 
and  linearly  polarized  electromagnetic  radiators,  where  the  far  field 
contains  only  a  single  component.  In  particular,  parallel  linearly 
polarized  current  elements  or  parallel  dipoles  are  covered. 

The  distant  field  in  an  arbitrary  direction  is  then  proportional 
to 

r  ib  +2«i(r  A)sin  0  cos(®-<a  ) 

V(0,*)  -  E(0,<p)  )  ak  e  k  k  ^  (1) 

k*l 


E(e,«p)  is  the  element  directivity  factor,  assumed  to  be  the  same  for 

all  elements.  The  phase  reference  has  been  taken  to  be  the  origin, 

0,  of  the  coordinate  system  and  the  factor  (rkA)  sin  0  cos^-^) 

is  simply  the  number  of  wavelengths  in  the  projection  of  0Pk  onto 

OD.  A  pencil  beam  with  a  maximum  in  the  direction  (0  ,®  )  will  be 

o  ,To 

formed  if  the  excitation  phases,  b^,  are  so  chosen  as  to  annihilate 
each  exponent  in  the  direction  (QQ,yQ).  This  choice  is 


bR  =  ^(r^X)  sin  0Q  cos(cpo  -  t^) 


(2) 


The  gain  in  any  direction,  (9,®),  is  the  ratio  of  intensity 
at  (0,®)  to  the  average  intensity  over  all  directions: 


0(0, tp) 


_ Iv(6,g)l2 

jb  |v(0,tp)  |2  dm 
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where  ...  dm  indicates  integration  over  the  sphere.  Gain,  G, 
without  qualification,  will  mean  G(0Q,cpo),  the  gain  in  the  direction 
of  the  beam  maximum.  This  gain  is  given  by 

|E(0o,<po)|2  |  E  ak|2/G  =  f  | sin  0|d8  f  |v(©,q>)  |2d*A* 
k=l  o  -n 


N-l  N  i(b  -b  )  N  p 

=  Re  E  E  2a  al  e  J  I  +  E  a^  L .  (3) 

j-1  k=j+l  J *  Jk  k=l  ^  kK 


where 


Xjk  "  J  ,Sin  0,t!I 

o  -n 


n  2rti 

e 


■|ln--^rk  cosfo-q^J-r  cos(®-<p  )] 


|e(0,(p)|  dq> 

00 


depends  only  upon  the  element  locations  and  the  element  directivity 
factor.  Notice  that 


I 


kk 


(5) 


is  merely  the  average  element  factor. 

When  the  element  factor,  E(0,cp),  has  the  twofold  synmetry 

|E(0,tp)|  =  |E(0,-p)|  =  |E(0,cp+n)|,  as  is  the  case  for  the  examples  to 

be  considered,  I  .  is  real  and  Eq.  (3)  becomes 
JK 


lE<0o’CTo>|2|  =«k1  N-l  N 


k=l 


=  2  E  S  a  a.  cos(b  -b  )  I 

J=1  k=J+l  J  K  k  j  JK 


N 

+  E 


kSL  ^  ^ 


(6) 
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SWAMPING  OF  ELEMENT  DIRECTIVITY  WHEN  THERE  IS  NO  THINNING 

The  gain  of  a  phased  array  will  change  as  the  beam  is  steered 
about.  Various  factors  contribute  to  this  change.  In  this  section 
evidence  is  given  that  the  element  factor  has  little  effect  on  the 
variation  of  gain  with  steering  whenever  the  array  is  producing  a 
sharp  beam  and  is  not  thinned.  In  Section  IV  these  assertions  are 
verified  by  comparing  the  variations  of  gain  with  steering  for  two 
arrays,  one  of  isotropic  elements  and  the  other  of  infinitesimal 
dipoles. 

The  quantity  |V(0,«p)|  of  Eq.  (3)  is  the  product  of  the  squared 
magnitudes  of  the  array  factor 

V  2«i(r  A)sin  0  cos({p-<a  ) 

A(0,tp)  =  ' ^  e 

k=l 

and  of  the  element  factor,  E(8,«p).  In  terms  of  A  and  E,  Eq.  (3)  may 
be  rewritten  as 

|E(0  ,<Pq)  A( 0q,«Po  )  1 2/G  =  f  I  sin  f  U(0,<p)  1 2 1  E(  0,«p)  |2d«p  (7) 

O  -  JT 

Suppose  the  element  directivity  is  a  slowly  varying  and  not  very 
peaked  function  of  position.  For  example,  if  the  elements  are 
parallel  infinitesimal  dipoles  the  element  factor  is  |E(0,<p)|  = 

sin  a,  where  a  is  the  angle  between  the  direction  (0,«p)  and  the  axis 
of  a  dipole.  This  pattern,  in  three  dimensions,  looks  like  a  smooth, 
holeless  doughnut.  When  the  pencil  beam  is  very  sharp  and  does  not 
come  near  a  null  of  the  slowly  varying  element  factor  and  when  most 
of  the  volume  Inside  the  array  factor  is  the  volume  directly  under 
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the  peak,  then  a  close  approximation  to  the  integral  in  Eq.  (7)  is 
obtained  by  setting  E(0,tp)  =  E( 0o,cpQ ) ,  its  value  at  the  peak  of 
A(9,«p),  whereupon  E(0o,(pQ)  cancels  out  of  Eq.  (7). 

Even  if  the  total  sidelobe  volume  of  A  is  not  negligible  there 
is  another  effect  that  reduces  the  dependence  of  gain  on  the  element 
factor:  the  contribution  to  the  integral  from  outside  the  peak  of 
A  stays  roughly  constant  as  the  beam  is  steered  around.  Denote  this 

2  2  p  2 

contribution  by  E^  A^(l  -  B),  where  E£  A^  is  an  average  value  outside 
the  peak  of  A  and  where  B  is  the  fraction  of  space  occupied  by  the 
peak  of  A.  Then  Eq.  (7)  may  be  written  as 


4  A*(l  -  B)  .  B  I*  A* 


or 


1 

G 


B  +  (1  -  B) 


-2  A2 
E1  *1 

2  2 

E  A 
o  o 


The  variation  of  G  is  produced  mainly  by  the  variation  of  Eq.  To 

the  extent  that  E£  A^/Eq  Aq  <  B  the  variation  in  the  second  term 

will  be  dampened  by  the  first  term,  B.  The  two  terms  will  be  about 

equal  when  half  the  energy  is  radiated  into  the  sidelobes,  but  even 

in  this  case  the  constant  B  will  roughly  halve  the  variation  of  G 

with  E  . 
o 

EFFECT  OF  THINNING  ON  GAIN  VARIATION  WITH  STEERING 

A  closely  spaced  array  approximates  a  continuously  excited  area. 
For  both  these  cases  gain  is  proportional  to  cos  0Q  for  angles,  0Q, 
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not  too  far  from  the  perpendicular.  For  widely  spaced  arrays, 

however,  this  cos  0Q  factor  will  disappear.  To  see  why  this  should 

be  so,  consider  the  two  parts  of  Eq.  (6).  1^  is  a  constant  and 

hence  the  second  term  is  independent  of  spacing  or  steering.  In  the 

first  term,  when  spacing  is  large  the  angles  b^  will  contain  many 

cycles  so  that  the  oscillating  factors,  cos(b  -  b  ),  are  equally 

k  j 

likely  to  be  positive  or  negative.  There  will  be  some  fluctuation, 
but  the  average  value  of  the  first  term  will  be  zero.  Thus  for  wide 
spacing,  gain  is  approximately 

N  .2 


iE<v*o>iX .  *  \ ) 


k-l 
N 

Ckk  E  “k 
**  k-l  K 


This  clearly  does  not  vary  as  cos  0 .  Since  I  ,  from  Eq.  (5),  is 

w  lUv 

the  average  of  |E(9,<p)|2,  it  follows  that 


l£(e  Q,tp Q)\2 

°  -VV’o) 


kk 


the  gain  of  the  element  factor  in  the  direction  (0Q,<po).  Thus  for 
wide  spacing  the  gain 


is  simply  a  constant  times  the  element  gain  in  the  direction  of 
steering. 
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The  argument  of  the  preceding  section  said,  that  element  gain 
should  be  swamped.  This  section  says  gain  follows  the  element  factor. 
However,  there  is  no  contradiction,  since  the  case  of  extreme  thinning 
is  precisely  the  case  for  which  E^  A^/E 2  A2  is  greater  than  B.  All 
these  remarks  are  checked  for  the  particular  examples  computed  in 
Section  IV.  In  summary,  when  spacing  is  close,  the  element  factor 
is  unimportant  and  gain  varies  with  steering  as  cos  0Q  out  to  a 
couple  beaawidths  away  from  the  equatorial  plane;  when  spacing  is 
great,  the  variation  of  gain  with  steering  follows  the  element  gain. 


E(e,<p)  = 


GAIN  IF  FT.FMENT  FACTOR  IS  UNIFORM  OVER  A  HEMISPHERE 
Let 

1  if  0  s  9  <  it/2 
0  if  n/2  <  e  <  n 

Such  a  radiator  will  be  referred  to  as  "isotropic"  for  short,  even 
though  its  effect  is  confined  to  a  half -space. 

From  Eq.  (4) 


rn//2  d9  pn 

=  !  sin  9  rr:  I 

o 


2ni  sin 


Sin’^-n 


■^"r^os^-t^J-rjCOsCtp-^j )] 


d (p 


,  M  r*  V 

.!  sin  0  ]-  .!  e 

o  -n 


dfp 


«/2 

=  f  sin  0  J  (2n(d  /X)  sin  0)  d0/2 
o  o  Jk 


(8) 
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sin(2«  <*1kA) 
XJk  4*  dJk/X 


(8) 

( cont'd) 


where 


A  '  rk  *  rj  '  2rjrk  cos(*k-*J)  ■  (vx/  *  (vy/ 


is  the  squared  distance  between  the  j-th  and  k-th  elements. 
When  J  =  k 


*kk  *  2 


(9) 


From  Eq.  (3)  the  gain  is  given  by 


r  N 

v  E  \  )  N-l  N  sin(2*d  /X)  .  .. 

—> —  A  k.A‘A  ~  *  ^  a  % 


do) 


An  important  special  case  is  ■  1;  then 


N-l  N 


N  i  i  *•-  *•  sin(2nd  /X) 

®  ‘ 2  * » A  k 


(n) 


Some  examples  for  this  case  appear  in  Section  IV.  For  wide  spacings 
the  terms  of  the  double  sum  go  down  because  of  the  divisor  d...  They 

JK 

also  tend  to  cancel  each  other  in  a  random  way  because  of  the  large 
and  irregular  values  of  (b. -b  ).  Hence  for  wide  spacings  G  <*>  2N. 

K  J 

The  factor  2  results  from  choosing  elements  that  radiate  "isotropically" 


into  a  hemisphere. 


12 


GAIN  FOR  DIPOLES  RADIATING  INTO  A  HEMISPHERE 

The  infinitesimal  dipoles  are  taken  to  lie  parallel  to  the  x- 

2  2 

axis.  The  element  factor  is  then  (l  -  cos  9  sin  0)  in  the  northern 
hemisphere,  0  s  9  s  «/2,  and  vanishes  in  the  southern.  For  the  calcu¬ 
lation  of  I,,  see  the  Appendix.  The  result  is 


"  sin  z 

-== —  -  cos  z 

^  sin  z  "I 

L  2  ' 

«.  Z 

^  z  J 

(12) 


where 


and 


cos  2w 


2  -  2*  a,kA 

K  - 


(13) 


(14) 


When  j  =  k 


I 


1 

kk  =  3 


(15) 


Gain  for  this  case  is  given  hy  Eq.  (3)>  vith  I ^  and  1^  given  by 

2  2  2 

Eqs.  (12)  through  (15 )  awi  with  |E(0Q,{po)|  =  1  -  sin  0Q  cos  <pQ. 

For  the  special  case  of  equal  excitation,  a]_  =  a2=  =  aN  =  ^"> 

the  gain  is  given  by 


13 


l-sin2e 


2 

cos  a 
o 

g7n 


1 

3 


2 

N 


N-l 

E 


H 

I 


j=l  k=j+l 


V  cos(VV 


(16) 


Some  examples  for  this  case  also  appear  in  Section  IV,  where 
the  behavior  of  gain  with  thinning  and  steering  is  compared  with  the 
corresponding  behavior  for  elements  with  uniform  radiation  over  a 


hemisphere. 


III.  DISCRETE  APPROXIMATION  OF  APERTURE  EXCITATIONS  BY  SPIRALS 


Consider  a  plane  array  intended  to  produce  a  steerable  pencil 
bean .  The  bean  shape  and  gain  will  depend  on  the  distribution  of 
elements  over  the  area  of  the  aperture  and  on  the  aperture  size.  The 
beaawidths  will  depend  mainly  on  the  aperture  diameters,  while  the 
relative  distribution  of  elements  will  determine  the  6idelobe  levels 
and  will  also  determine  whether  spurious  main  lobes  arise.  The 
general  behavior  of  near  sidelobes  can  be  controlled  by  making  the 
element  distribution  approximate  a  continuous  excitation  having  a 
known,  and  acceptable,  set  of  sidelobes.  For  example,  a  suitable 
continuous  excitation  can  be  divided  into  areas  of  equal  total  ex¬ 
citation  and  an  element  of  a  discrete  array  assigned  to  each  such 
area. 

The  actual  partitioning  of  the  illumination  into  equal  incre¬ 
ments  can  be  done  in  many  ways.  For  example,  a  square  aperture  can 
be  divided  into  a  regular  array  of  rectangles  and  an  element  assigned 
to  each  rectangle,  if  a  uniform  excitation  is  to  be  approximated. 

When  such  a  regular  array  is  used,  spurious  lobes  arise  as  soon  as 
the  spacing  exceeds  one  wavelength.  By  randomizing  the  spacings, 
such  lobes  can  be  decreased.  The  general  objective  is  so  to  space 
the  elements  that  the  distribution  of  their  projections  on  any  line 
is  always  irregular.  This  immediately  suggests  basing  the  distri¬ 
bution  on  a  set  of  uniform  distributions  on  circles,  since  circles 
already  are  free  from  preferred  directions.  If  an  annular  region  is 
to  be  occupied  by  a  fixed  number  of  elements,  it  would  presumably  be 
divided  evenly  and  an  element  assigned  to  each  equal  part,  pe.haps 
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being  randomly  displaced  inside  its  assigned  area.  The  rings  could 
be  rotated  with  respect  to  one  another  so  that  no  strongly  lined-up 
elements  formed  spokes. 

From  the  above  discussion  it  is  only  a  short  jump  to  the  idea  of 
dividing  the  continuous  excitation  into  equal  pieces  by  winding  out 
from  the  center  in  a  spiral  strip.  The  spiral  strip  will  then  be 
divided  into  roughly  square  pieces,  each  to  be  represented  by  one 
element  of  the  final  array.  If  the  array  is  to  approximate  a  con¬ 
stant  excitation,  then  the  spiral  strip  will  soon  reach  a  nearly 
constant  width.  If  the  illumination  is  to  be  tapered,  then  the  strip 
should  gradually  widen  as  the  spiral  expands.  Similarly,  the  spacing 
of  elements  along  the  strip  would  be  even  if  a  constant  illumination 
is  to  be  approximated,  but  would  gradually  increase  if  there  is  to 
be  a  tapering  at  the  edge. 

There  is  no  guarantee  that  such  a  spiral  plan  for  positioning 
the  elements  will  not  produce  spokes,  or  even  nearly  parallel  lines 
of  elements.  However,  it  seems  probable  that  such  accidents  would 
be  rare. 

DISC  WITH  CONSTANT  EXCITATION 

Consider  the  problem  of  piewing  a  given  number  of  elements 
inside  a  circle  so  that  the  spacing  between  elements  is  roughly 
uniform  but  so  that  the  elements  do  not  tend  to  fall  on  straight 
lines.  Suppose  the  elements  are  evenly  spaced  along  the  arc  of  an 
expanding  spiral.  After  the  spiral  has  grown  a  bit,  the  distance 
between  successive  elements  on  the  spiral  (chord  or  secant)  will 
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soon  be  very  nearly  the  sane  as  the  distance  along  the  arc,  hence 
choosing  equal  spacing  along  the  arc  guarantees  roughly  uniform  spacing 
of  elements  in  the  9  direction.  The  direction  of  the  discrepancy  will 
be  such  as  to  produce  a  slight  tapering  down  as  the  outside  of  the 
array  is  approached,  which  is  a  discrepancy  in  the  right,  or  helpful, 
direction  (since  it  means  slightly  lower  sidelobes).  To  determine 
the  spiral  approximately,  consider  the  relation  between  arc  length, 
s,  and  radius,  r.  Let  d  be  the  constant  distance  between  elements 
along  the  arc  of  the  spiral;  then  the  number  of  elements  in  length 
ds  will  be  ds/d.  Suppose  the  radius  grows  by  dr  during  the  space  ds. 
If  a  constant  illumination  is  to  be  approximated,  then  the  number  of 
elements  in  an  annulus  of  width  dr  should  be  proportional  to  the  area, 
2:trdr,  of  the  ring.  Hence  the  spiral  is  given  by 

ds  =  r  dr/rQ  ^7) 


where  r  is  a  proportionality  constant.  Integrating  Eq.  (17)  results 
o 

in 


p  2 

r  a  r  +  2r  s 
o  o 


(IS) 


From  Eq.  (IB)  and  ds2  -  r2  d02  +  dr2  it  follows  that 


d0 


r  +  2r  s 
o  o 


or 


(19) 
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Equations  (l8)  and  (19)  are  parametric  equations  of  the  spiral. 
Spacing  elements  equally  along  the  arc  means  that  the  k-th  element 
corresponds  to  parameter  sk  =  (k  -  l)d,  for  k  =  1,  2,  ...,  N,  with 
the  polar  coordinates  of  the  k-th  element  given  by  Eqs.  (lB)  and 
(19).  Let  the  outermost,  or  N-th,  element  lie  at  distance  a  from 
the  center;  thus  a  may  be  considered  the  radius  of  the  roughly 
circular  array.  Then 

a2  -  r^  +  2rQ(N  -  l)d  (20) 

A  choice  for  the  starting  radius,  rQ,  has  still  to  be  cade.  The 
area  represented  by  an  individual  element  will  be  roughly  square  if 
the  distance  between  successive  turns  of  the  spiral  is  equal  to  the 
distance,  d,  between  elements  along  the  spiral.  This  will  make  the 
radial  and  tangential  components  of  the  interelement  distances  about 
the  same.  An  approximate  solution  for  this  requirement  is  possible 
if  the  spiral  has  several  turns.  Near  the  outside  the  term 

ySs/i-Q  is  approximately  n/2  and  can  be  ignored.  From  Eqs.  (lB) 
and  (19),  for  two  values  of  the  parameter,  t  and  s,  corresponding 
to  a  change  of  2«  in  0 
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hence 


and 


d  as  2«  rQ 


(21) 


The  elements  are  located  where  sk  =  (k  -  l)d.  Take  the  N-th 
element  at  the  outer  radius,  a;  then 

a2  =  r2  +  2r  (N  -  l)  2k  r 
o  o  o 


or 

r  -  ,  a  (22) 

0  's/I  +  4k(N  -  1) 

EXAMPLE- -100  ELEMENTS  OVER  A  DISC 

If  the  outer  radius  is  a  =  l4.1,  then  rQ  =  l4.l/ 1+4k  ( 100-1 ) 

„  .400  and  d  =  2«  rQ  =  2.51.  From  Eqs.  (lB)  and  (19)  there  result 

rk  2.000k  -  1.848  (  23) 

and 

0k  =^J  12-57(k  -  1)  -  tan"\/l2.57U  -  l)  (2»0 

The  100  locations  for  this  example  are  plotted  in  Fig.  2.  They 
look  uniformly  distributed  and  not  lined  up  in  any  particular 


direction. 
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Spiral  distribution  approximating  uniform  excitation 
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It  remains  to  discover  the  extent  to  which  this  arrangement  is 
"random"  and  has  succeeded  in  avoiding  big  si delobe s  when  the  average 
interelement  spacing  is  made  noticeably  greater  than  half  a  wavelength. 
The  pattern  produced  by  100  elements  arranged  as  in  Eqs.  (23)  and 
(24)  is  described  in  Section  IV,  along  with  the  other  numerical  re¬ 
sults  for  this  example. 


IV.  A  DETAILED  ILLUSTRATION 


To  verify  the  conclusions  of  Section  II  an  example  is  worked  out 

in  detail.  The  variation  of  gain  with  steering  and  spreading  (or 

thinning)  is  found  for  the  spiral  array  of  Fig.  2,  using  both  isotropic 

* 

elements  and  infinitesimal  dipoles.  For  no  spreading,  the  gain 
varies  with  steering  according  to  the  well-known  cos  6q  approximation, 
while  for  great  spreading  the  variation  of  array  gain  with  steering 
follows  the  gain  of  the  element  factor. 

THE  THINNING  PARAMETER,  S 

Up  to  now,  thinning  or  spreading  has  been  discussed  in  a  rough 
quantitative  way  only.  Before  any  calculation  can  be  done,  the  notion 
of  thinning  needs  to  be  made  exact.  For  a  regular  array  whose  elements 
lie  at  the  intersections  in  a  square  grid,  the  spacing  parameter  would 
naturally  be  defined  to  be  the  mesh  dimension  in,  perhaps,  half  wave¬ 
lengths.  For  the  sort  of  irregular  array  presented  by  the  spiral  in 
this  example,  no  such  direct  and  simple  definition  is  possible. 

Ary  definition  of  the  spacing  parameter,  S,  will  always  be 
somewhat  arbitrary.  The  one  to  be  given  here  has  two  merits:  it  is 
fairly  natural,  and  it  reduces  to  the  "right"  number  for  regular, 
rectangular  arrays.  Let  A  be  the  area  of  an  array  which  is  occupied 

- 

Because  there  wan  assumed  to  be  no  radiation  into  the  southern 
hemisphere,  the  results  may  be  3  db  higher  than  might  be  expected. 

No  special  virtue  inheres  in  infinitesimal  dipoles;  they  just  happen 
to  be  a  convenient  example  for  which  to  verify  the  assertions 
about  the  effect  of  the  element  factor  on  the  variation  of  gain  with 
thinning  and  steering. 
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by  N  elements.  Then  the  spacing  (or  thinning,  or  spreading) 
parameter,  S,  is  defined  as 

S  -  — ^  W 

H(A/2 y 

The  units  are  square  halfwaves  per  element.  A  regular  array  of  N 
elements,  spaced  K  halfwaves,  has  S  ma  K^,  with  =  S  =  1  the  normal 
or  unthinned  case.  For  an  array  which  has  been  literally  thinned 
(by  removing  elements  more  or  less  uniformly,  while  leaving  the  same 
area  occupied),  S  is  the  factor  by  which  the  number  of  elements  in 
the  array  has  been  reduced  from  the  normal  case  of  halfwave  spacing. 

A  question  still  to  be  settled  is  what  area  does  the  spiral  of 
Fig.  2  actually  "occupy."  For  the  results  to  be  exhibited,  the  area 
occupied  is  taken  to  bn  a  circle  whose  diameter  is  half  a  wavelength 
greater  than  the  diameter  of  the  spiral.  The  diameter  of  the  spiral 
is  the  maximum  distance,  D^,  between  a  pair  of  elements  on  the 
spiral.  Adding  the  extra  half  wavelength  reflects  a  belief  that  the 
sphere  of  influence  of  a  simple  radiator  like  a  dipole  extends  out¬ 
ward  about  a  quarter  wavelength.  The  spreading  factor  for  the  spiral 

array  is  then 

*<V.  * 

s  =  - ==-5 - 

N  A 

VARIATION  OF  GAIN  WITH  THINNING 

Figures  3  and  4  are  almost  self-explanatory.  For  a  beam  steered 
from  the  north  pole  to  the  equator  along  longitude  0°  the  gain 
(normalized  by  N,  the  number  of  elements)  is  plotted  vs  0  for  values 
of  the  thinning  parameter,  S  «  !,•••,  8*  Figure  3  is  for  so  called 


S  *  square  halfwaves  per  element 
-  0 

Fig.  4  — 


versus  6° 


iir 
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i60tropic  elements  (uniform  radiation  over  hemisphere),  while  Fig.  4 
is  for  dipoles  except  that,  again,  radiation  into  the  southern 
hemisphere  was  assumed  zero. 

To  show  how  well  the  cos  6  law  is  followed  for  the  case  S  =  1, 
the  3  db  limit  of  the  endfire  bean  has  been  sketched.  Since  the 
cos  0  rule  would  predict  zero  gain  at  0  =  90°,  it  obviously  cannot 
hold  over  the  whole  range.  The  example  shows  that  the  rule  is  closely 
followed  out  to  within  a  beanwidth  (endfire)  of  90°. 

The  approximation,  cos  0,  is  shown  for  the  unthinned  case 
(S  =  1),  and  the  element  factor  is  shown  for  the  most  thinned  case 
(S  =  8).  For  dipoles  along  OX,  the  element  gain  varies  as  cos2  0, 
which  is  the  element  factor  plotted  in  Fig.  4  for  S  =  8. 

Figure  5  contains  the  same  data  as  Fig.  3,  but  plotted  against 
spacing  factor  for  constant  steering  angle,  0. 

Figure  6  exhibits  the  variation  of  gain  for  a  wide  range  of 
s pacings  when  the  beam  is  aimed  broadside  (0=0°).  The  solid  curve 
is  for  isotropic  elements  while  the  circles  are  for  infinitesimal 
dipoles.  Note  that  at  low  S  the  gains  are  virtually  identical,  while 
as  S  increases  the  gains  gradually  come  to  differ  by  about  1.8  db, 
corresponding  to  the  gain  (3/2)  of  a  dipole. 


COMPARISON  OF  SPREADING  AND  TWTTJNING 

All  the  results  plotted  in  Figs.  3  through  6  have  been  for  a 
fixed  number  (100)  of  elements.  They  therefore  show  the  effects  of 
spreading  a  fixed  number  of  elements  over  a  greater  and  greater  area 
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keeping  the  configuration  similar  to  the  original.  What  if  the  area 
were  to  be  held  fixed  and  a  given  distribution  actually  thinned  by, 
say,  halving  and  quartering  the  number  of  elements?  The  answer  to 
this  question  for  the  spiral  example  appears  in  Fig.  7-  The  three 
curves  on  the  left  are  repeated  from  Fig.  3-  The  three  curves  on 
the  right  also  correspond  to  spacing  factors  of  S  =  1,  2  and  4,  but 
these  factors  were  produced  by  reducing  the  number  of  elements  rather 
than  by  spreading  100  elements  over  a  greater  area. 

The  cases  considered  fall  in  the  most  critical  region  of  spacing- 
the  transition  region  between  no  thinning  and  great  thinning.  Never¬ 
theless,  the  variations  of  gain  with  steering  for  these  two  kinds  of 
thinning  (increasing  the  area  vs  decreasing  the  number  of  elements) 
are  in  remarkable  accord. 

FPCTTnQM  OF  ARRAY  FACTOR  FROM  GRATING  LOBES 
The  detailed  array  factor, 

A(0,tp)  =  E  exp((2«iA)rk[sin  0  cos(«p-<Pk)  "  3*n  eQ  c0S(V*k^ 
k=l 

for  an  array  of  equally  excited  elements  may  be  written  in  either  of 
the  universal  forms 

A(U,V)  =  E  exp((2«i/X)((xk//3)U  +  (yk/Vs)vl} 

k=l 


or 


A(Z,Y)  =  2  expf (2rtiA)(rk//S)Z  cosCY-qj^ 

k=l 


100  Elements 
variable  area 


Area  *  100  square  halfwaves 
variable  number  of  elements 


Fig.  7  — Variation  of  gain  with  steering  and  thinning 
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where 

U  -  Z  COS  T  -  /5(sln  «  cos  »  -  sin  90  cos  *0) 

V  =  Z  sin  I  =  /§(sln  9  sin  9  -  sin  6Q  sin  mQ) 

The  scsling  xJ/3,  y^/S,  r^/3  pioduces  n  normalised  set  of 
element  locations  for  which  the  new  spacing  Is  one  square  halfcave 
per  element.  A  universal  array  factor  rmy  be  computed  In  terms  of 
the  pair  of  variables  (U,V)  or  the  pair  <M).  These  pair,  esibody 
both  steering  and  spreading  variations.  If  3  -  1  and  the  barns  Is 

aimed  brtmdsid.  (.„  -  0),  then  |z|  «  1  «*  -1  ‘  4  1 

possible  points  In  real  space.  If  there  Is  to  be  steering  do»n  to 
,  .  s/2,  then  the  «nges  of  «,V  and  Z  must  be  doubled  to  aecount  for 

all  of  reel  space  and  for  steering  as  veil.  Finally,  If  there  1.  to 
be  thinning  by  the  factor,  S,  then  the  ranges  of  U,V  and  Z  are 
multiplied  by  a  further  /%.  For  exmsple,  If  the  array  of  the  example 
la  to  be  fully  steered  and  If  thinning  by  a  factor  6  .  16  is  to  be 
used,  then  the  range,  of  U,V  .nd  Z  are  0  s  Z2  ■  U2  *  V2  s  6*  or 
0  s  z  s  8.  Several  slice.  (I  -  0°,  30°,  60°,  90°)  »cre  plotted  for 
the  array  factor  of  the  example  and  are  shorn  In  Figs.  6  mid  9-  In 
addition  a  coarse  survey  of  the  array  factor  vas  male  by  computing 
It  at  a  101  by  SOI  grid  of  points  In  the  reetsngular  Interval 
0  s  U  s  50.  -50  s  V  s  50.  The  highest  value  encountered  va.  10.7  db 
dovn  from  the  main  beam.  Allowing  a  factor  of  S  for  steering.  It 
appears  that  this  100-element  configuration  may  be  free  of  grating 
lobes  even  when  average  Interelesmnt  spmrlng  Is  IS  wavelengths,  though 
dr,,!  .nm.  were  not  dense  enough  to  guarantee  that  some 


the  survey 


i 


l 


Fig. 9__ slices  of  space  factor  of  100-element  spiral 
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lobe  has  not  been  missed.  The  survey  amounts  to  t-ain ng  about  one 
point  for  every  three  or  four  lobes,  the  samples  being  equally  spaced 
while  the  lobe  s pacings  are  randan. 

The  purpose  of  this  discussion  is  not  to  suggest  that  the  par¬ 
ticular  example  presented  is  an  outstanding  array  (after  all,  it  has 
17  db  nearby  sidelobes),  but  rather  to  illustrate  the  success  of  the 
spiral  layout  in  avoiding  grating  lobes  when  spread  or  thinned  to  an 
average  interelement  space  of  2 X,  with  full  steering.  As  a  comparison, 
notice  that  a  10  by  10  square  array  with  full  steering  will  produce 
a  grating  lobe  at  X/2  spacing  (at  endfire  it  is  bi-directional),  and 
that  its  first  sidelobe  is  12  db  down.  The  17  db  first  sidelobe  of 
the  example  is  about  what  one  expects  from  o  uniformly  excited 


circular  area. 


& 
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Appendix 


FOR  INFINITESIMAL  DIPOLES 


2 

The  element  factor  for  an  Infinitesimal  dipole  is  sin  A  (in¬ 
tensity),  where  A  is  the  angle  between  the  dipole  axis  and  a  vector 

to  the  distant  point.  If  the  dipoles  are  all  parallel  to  the  x-axis, 

2  2  2 
then  sin  A  -  1  -  sin  9  C06  tp  and 


I  -  r|8in0|d9J  e  X  L  k  J  J  fl-sin  Ocos  ojfc 

A  "  ^ 


The  first  step  in  evaluating  I  is  to  rewrite  the  bracket  in  the 
exponent  as 

rkcos(g-^)-rjCos(tp-<pj )  = 


where 


Jk 


■  i  * r]  ■  *3  v05 


(cV^} 


(VX/  +  (Vy/ 


rkco5^  •  r1COGC?j 

006  v - 


rk6inc^  -  r.isiIV?j 

ein  *.k  =  - - 

Jk 


Shifting  the  c  -  integral  by'  an  amount  gives 
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n 

2 


R  2kL  sine  . 

IJk  =  2  J  sinede  J  e  x  Jk 


o  -it 

it 


it  2ni  sine 


a  J  sin^ede  j  e-*--  V0**^ 


-It 


(9+*jk)d<P 


Now  put  2coe2(9+tjk).  1  +  cos2cp-cos2*Jk  -  sinapBin2ifrjk  and  split  one 
of  the  sin3e  factors  into  sino(l-cos2e).  The  result  is 


XJk-a 


n 

2  n  2«1  sing  ,, 

J  sinede  J  e  x  ^dcp 


-it 


n 

2 


n  2}ii  sine 


-  J  (sine-sine  cos2e)de  J  e  x  jkCOE<? 


dtp 


-it 


2  k  2d  sine  , 

~ —  Vo*p 


-  cos2^k  J  Bin3ede  J  e 


cos2cpd<p 


-R 


2  r  Bine  , 

Sin2*Jk  |  Bln39de  J  «  X  JkC0SCPBin^do 


The  last  integral  vanishes  because  the  integrand  is  an  odd  function 
of  9.  The  first  two  integrals  simplify  a  bit,  so  that  finally 


I 


* 

2 


Ijk  =  2*  J  sine  J0(2«djjt8ineA)de 


« 

2 


+  2*  |  sine  cog  9  J0(2«djk6ine/x)de 


« 

2 


2jtcos2*  J  sin3e  J2(2jtdJksine/\) 


Use  has  been  made  of  the  fact  that 


i  n  - 

i  |  eizco6*  cos  ncp  dtp  =  e  °  2  Jjz) 
o 

The  three  integrals  in  IJk  are  special  cases  (see  Ref.  11,  p.  373)  of 
B(a,b,z)  =  J  Ja(zsine)sinB+10  cos^^e  do  =  2  £ii+1"  Ja+b+l(z) 


Particular  cases  needed  are 


a  =  o  ,  o  =  2 


z  cos  z 


a  =  o  ,  b  =  -  5 


1  \  r(l/2)  T  _  sin  z 

B(o>-  2>Z)  =  TTT  J1/2(Z)  =  ~~ 


I 
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and 


a  -  2,  b  =  -  ± 


B<2'-  i**>  ■  ^  J5/2w  J\  ( 


sin  z  _ N  sin  2 

- r -  “  COS  2  , 


Let  zJk  ■  2«djk/X;  then 


V  cln  fjk  , ,la>  c°*,> 


Jk 


Jk 


Jk 


cos2*jk  [  4-  C 


sin  2 


Jk 


sin  2 


cos  z 


ZJk  "  Jk 


Jk 


Jk 


This  expression  is  indeterminate  f or  J  =  k  (and  hence  z  =  o);  how- 

Jk 

ever,  its  limit  as  z  .  o  is 

JK 

T  -  8x 

*kk  -  T 


This  value  is  also  obtained  directly  from  the  basic  formula: 

it 

2  « 

Ij^  =  2  j  Sine  de  |  (1-Sin2e  cos2cp)d<p  =  ^ 
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